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The dependence of the electrical conductivity on the oxygen partial pressure is calculated for the prototypical perovskite 
BaTi03 based on data obtained from first-principles calculations within density functional theory. The equilibrium 
point defect concentrations are obtained via a self-consistent determination of the electron chemical potential. This 
allows to derive charge carrier concentrations for a given temperature and chemical environment and eventually the 
electrial conductivity. The calculations are in excellent agreement with experimental data if an accidental acceptor 
dopant level of 10 17 cm -3 is assumed. It is shown that doubly charged oxygen vacancies are accountable for the high- 
temperature n-type conduction under oxygen-poor conditions. The high-temperature p-type conduction observed at 
large oxygen pressures is due to barium vacancies and titanium-oxygen di-vacancies under Ti and Ba-rich conditions, 
respectively. Finally, the connection between the present approach and the mass-action law approach to point defect 
thermodynamics is discussed. 
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I. INTRODUCTION 

Point defects control the functional properties of semicon- 
ductors and many insulators, but are usually difficult to as- 
sess experimentally. If bulk properties like conductivity or 
diffusivity are measured, it is necessary to introduce model as- 
sumptions in order to relate the measured macroscopic quan- 
tities to microscopic point defect properties. In contrast, lo- 
cal probes such as electron spin resonance or positron anni- 
hilation spectroscopy can provide very specific information 
on point defect structures but are usually restricted to certain 
electronic configurations (unpaired spins) or types of defects 
(open volumes, vacancies). In general, only by combining 
several experimental probes a consistent description of the de- 
fect chemistry in a given material can be obtained^ To compli- 
cate things further, the correlation between defect properties 
and any experimentally measured response is typically indi- 
rect and often prone to ambiguities. 

Computational modeling techniques evolved rapidly in re- 
cent years, in particular in the realm of first-principles cal- 
culations. Among these methods schemes based on density- 
functional theory (DFT) are extremely popular as they have 
become increasingly reliable and have been shown to be capa- 
ble of predicting various materials properties. First-principles 
modeling is particularly attractive with regard to point defects. 
It allows to obtain detailed information about thermodynamic 
and kinetic properties (formation energies and volumes, mi- 
gration barriers and entropies) as well as the electronic struc- 
ture which — at this level of detail — are not available through 
experimental techniques (see e.g., Refs. However, only 
in very few cases calculations of point defect properties have 
been employed to derive macroscopically measurable quanti- 
ties such as conductivities or diffusivities (see e.g., Refs. and 
7). It is however instrumental to develop these connections 
between calculation and experiment in order to verify the un- 
derlying methods and to establish their predictive power. 



In the present contribution we demonstrate for the case of 
BaTi03 how theoretical data obtained from first-principles 
calculations 5 can be used to derive the dependence of the elec- 
trical conductivity on the oxygen partial pressure. The electri- 
cal conductivity is a technologically highly relevant property, 
the understanding of which is at the very foundation of de- 
vice technology. The successful modeling of this property de- 
scribed in this work provides the basis for future work which 
should address, e.g., the role of kinetic effects, extrinsic de- 
fects, or defect association. 

Barium titanate is a prototypical ferroelectric material with 
a paraelectric-ferroelectric transition temperature of 393 K. 
Its most important technological application is in thin-film 
capacitors^ In addition BaTiOs serves as an end member in 
several lead-free ferroelectric alloys,— and is used — often in 
combination with SrTiOs — to obtain tunable RF devices 
Because of its technological importance it has been exten- 
sively investigated, and a reliable as well as extensive data 
base is available (see e.g., Refs. [T2I - I20I) . BaTiOs is therefore 
not only a very interesting material for theoretical investiga- 
tions, but it provides also an excellent testbed for carrying out 
a stringent comparison between calculated and experimental 
data. 

In the following we first introduce the thermodynamic 
framework and the relevant equations of semiconductor 
physics. Combining these equations we are able to determine 
self-consistently for a given temperature and a given chem- 
ical environment (1) the electron chemical potential, (2) the 
point defect concentrations, (3) the charge carrier concentra- 
tions, and eventually (4) the electrical conductivity. In section 
Sect. |nl] the results of the calculations are compared to high- 
temperature experimental data and the dependence of the con- 
ductivity on the oxygen partial pressure is analyzed in terms 
of the point defect equilibria in the material. The model is 
subsequently employed to extrapolate the materials behavior 
to lower temperatures, where experimental measurements are 
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no longer available. 

II. THERMODYNAMIC FORMALISM 

A. Gibbs free energy of point defect formation 

The Gibbs free energy of formation of a point defect in 
charge state q can be consistently derived from thermody- 
namic principles and depends on the chemical potentials /Ltj 
of the constituents ("the chemical environment") and the elec- 
tron chemical potential [i e as follow o 21 ' 22 

AG d = (G def - G host ) - Y Ar W + l( E VBM + Ve), 

j 

(1) 

where Gdef and Ghost are the Gibbs free energies of the sys- 
tem with and without the defect, respectively. The difference 
in the number of atoms of type i between these two systems 
is denoted Anj and the sum runs over the elements present in 
the system. It is convenient to separate the chemical poten- 
tial into the chemical potential of the ground state fj^ and the 
variation relative to the ground state chemical potential A/Zj, 
i.e. 

Mj^Mj+A/ij. (2) 

Finally, the position of the valence band, £Vbm, defines the 
reference of the energy scale for the electron chemical poten- 
tial, /i e . 

Knowledge of the defect formation energy allows to cal- 
culate the defect concentration as a function of temperature, 
which is given by 

where c° denotes the number of sites available for defects on 
the respective sublattice per volume (e.g., the density of bar- 
ium sites in the case of barium vacancies). 

B. Charge neutrality condition: the intrinsic electron 
chemical potential 

The electron chemical potential fj, e , which appears in 
Eq. (0, is actually not a free parameter but fixed by the charge 
neutrality condition 1 , 

n e + n A = n h + n D , (4) 

which links the concentration of intrinsic electrons n e and 
holes rih to the concentration of charge carriers induced by 
acceptors ua and donors ud- As discussed in the following 
each term in Eq. (0]i is exponentially dependent on the electron 
chemical potential fi e . Finding a solution of Eq. © therefore 
yields the intrinsic (self-consistent) chemical potential for a 
given temperature and chemical environment (Sect. HTDl . 



The intrinsic charge carrier concentrations are obtained by 
integrating the number of unoccupied states up to the valence 
band maximum (VBM) and the number of occupied states 
above the conduction band minimum (CBM), 

n e (/i e ) = / D(E)f(E,n e )dE (5a) 

JCBM 

/VBM 
D(E)[l-f(E,n e )]dE, (5b) 
-OO 

where f(E,n e ) = {1 + cxp [(E - fieJ/taT]}^ 1 is the 
Fermi-Dirac distribution. 

The concentrations of point defect induced carriers are 
obtained by summing the concentrations of acceptors and 
donors, 

donors , ^ d v 

no = Y e * c ° ex P ( ) (6a) 

i \ B / 

acceptors . ^ d . 

i \ B J 

where e is the unit of charge and qi is the charge state of de- 
fect i. Additional charge carriers contributed by dopants or 
impurities ("accidental dopants") can be simply added to nA 
and no, respectively. Their concentrations are given by 

ng l = eq D c^ ext [1 -f(E G -E D , Me )] (7) 

nf = eq A c A ext [f(E A ,f, e )}, (8) 

where Ed and E A are the donor and acceptor equilibrium 
transition levels measured with respect to the conduction band 
minimum (CBM) and the valence band maximum (VBM), re- 
spectively, and c° D ,ext and c° A ,ext are the impurity concentra- 
tions. 

Mathematically, it is possible that there is more than one 
solution of the charge neutrality condition Eq. This can 
occur e.g., if some defect concentrations change as a func- 
tion of temperature whereas some others are held constant. In 
such a case the solution with the lowest Gibbs free energy is 
selected. The difference of the Gibbs free energy with respect 
to the equivalent defect-free reference system is obtained by 
summing the formation energies of all defects in the system 
minus the configurational entropy, 

&G!*Y2 NiAG i~ kBT]nSl > (9) 

i 

where denotes the number of defects of type i and £1 de- 
notes the number of possible configurations (compare chap- 
ter 3.3 of Ref.H). 

C. Electrical conductivity 

The electrical conductivity is obtained by summing over all 
mobile charge carrying species^ 

a = Y B i c Mi e ( 10 ) 
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where Bi are the mobilities, c, are the concentrations per vol- 
ume, and qi are the charge numbers. Typically one distin- 
guishes the electronic and ionic contributions a = a e i + <Ti on . 
According to Eq. (TTOb the former is simply 



a e l = BeTLee + BhJlhe 
elections holes 



(ID 



where n e and are given by Eqs. Q under the constraint 
that the charge neutrality condition is fulfilled. Charge car- 
rier mobilities subsume the contributions of all possible scat- 
tering mechanisms — most importantly defects and phonons — 
and are therefore very difficult to calculate. At present we re- 
sort to experimental data instead. For the electron mobility 
we use the expression given in Ref. 14 which is a fit to single 
crystal data from Seuter using the expression given by Ihrig 2 ^ 



orn 2 K 3/2 

Be = 8080 • T- 3 / 2 • exp 



0.021 eV 



k B T 



(12) 



For the hole mobility we follow Ref. |l_4j and assume Bh « 
B e /2. 

In order to obtain the ionic conductivity one can use 
the Einstein-Smoluchowski relation Bi = eDi/ksT to re- 
place the defect mobility with the defect diffusivity Di = 
D exp (— AG™/fc s T) which yields 



Eq i e 2 D c l 
— r m ex P 



AG" 1 



(13) 



For cubic crystals the pre-factor is D = 6r aQ where ao is 
the lattice constant and Tq is the attempt frequency. The latter 
can be approximated by the lowest optical phonon frequency 
which yields 24 To « 5THz and D$ w 10~ 3 cm/s 2 . The mi- 
gration energies for intrinsic vacancies have been reported in 
Ref. |H Since both the migration entropy and the migration 
volume are about a factor of magnitude smaller than the for- 
mation entropy and volume, they can be safely neglected in 
the present case, i.e., we can assume AG™ « AE™ ■ Using 
these data it is found that for the present material the ionic 
contribution at elevated temperatures is about four orders of 
magnitude smaller than the electronic contribution. In the fol- 
lowing we therefore consider the electronic part only. 



D. Phase stability: limitations on the chemical potentials 

The chemical potentials, = + A/ij, which appear 
in Eq. ([T]), are subject to several thermodynamic constraints. 
First, they cannot become more positive than the chemical po- 
tential of the reference phase, i.e. Afij < 0, where the ref- 
erence phase for oxygen is the 02 molecule, for barium the 
body-centered cubic crystal, and for titanium the hexagonal- 
close packed crystal. If any chemical potential reaches its 
upper limit, the respective elemental ground state phase pre- 
cipitates. Second, the chemical potentials of the constituting 
elements are coupled by the requirement that 5 

A/iBa + A/ixi + 3A^ = Ai?/[BaTi03], (14) 
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FIG. 1 . Phase diagram for cubic barium titanate at zero Kelvin as 
determined from density functional theory calculations (Ref.Q). The 
area confined between points A, B, C and D is the chemical stability 
range of BaTi03. The chemical potentials confined to the lines A-D 
and B-C are referred to as B a and Ti-rich in the text, respectively. 



where AiJ/[BaTi03] is the formation energy of BaTi03. 
Further constraints result from the formation of competing 
phases, namely 



A^Ba + A/i < A#/[Ba0] 
A/ij; + 2 A/i0 < AH f [Ti0 2 ] 



(15a) 
(15b) 



If all of these restrictions are included, one obtains the static 
phase diagram for T=0 K depicted in Fig. Q] The outer tri- 
angle follows from condition ( TBI while the lines separating 
the BaTiC<3, TiC<2 and BaO phases result from Eqs. (1151 . The 
gray shaded area is the (zero Kelvin) stability range of BaTiOs 
with respect to BaO and TiC<2. 

Experimentally, the way to control the thermodynamic 
boundary conditions is to use either BaO or Ti02 excess dur- 
ing materials processing, and to vary the oxygen partial pres- 
sure, P0,, during processing and measurements. Adjusting 
the excess of either Ba or Ti corresponds to constraining the 
accessible range of chemical potentials to the lines A-D (Ba- 
rich limit, equilibrium between BaO and BaTiOs) or B-C (Ti- 
rich limit, equilibrium between Ti02 and BaTiOs) in Fig. Q] 
Varying the oxygen partial pressure (i.e., to the oxygen chem- 
ical potential) is equivalent to moving along these lines where 
the extremal points A and B on one, and C and D on the other 
side correspond to metal-rich (low p©,) and oxygen-rich (high 
P0 2 ) conditions, respectively. 

While the formation energies calculated in Ref. HI are given 
as a function of the chemical potentials, experimentally the 
conductivity is measured as a function of the oxygen partial 
pressure. In order to compare the conductivities as calculated 
for different chemical potentials with experimental data, one 
must therefore convert between the oxygen chemical potential 
and the oxygen partial pressure. The two quantities are related 
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according ti 



,25.26 



V0(T,P0 2 ) = ^ (T,p Q ) + -k B T\n 



Pok 



^ (16) 



where p^ 2 denotes the reference pressure. We choose the 
isolated oxygen dimer molecule as the zero Kelvin reference 
state, jLt0(OK,j5°) = \Eq 2 . For consistency with the exper- 
imental data and following Ref. |25| we use the experimental 
value for E<$ 2 = — 5.16eV/dimer and the experimentally de- 
termined temperature dependence of /i0(X \P0 ) (Ref.l27l). 

The phase diagram in Fig. Q]is strictly valid only at zero 
temperature. At finite temperatures the construction would 
have to be based on the free energies of formation instead. 
The major effect arises from the differences between the vi- 
brational entropies between the various relevant phases. With 
the exception of oxygen all these phases are crystalline. The 
entropies of crystalline solids are, however, much smaller than 
the entropies of gases. By the far the most important term is 
therefore the change of the free energy of the oxygen reser- 
voir, which is properly taken into account via Eq. ( TToT l. It is 
therefore admissible to use the phase diagram established here 
also at finite temperatures. 



E. Summary of algorithm 

In summary computing the conductivity proceeds as fol- 
lows: (/) The electron chemical potential is self-consistently 
determined as described in Sect. Ill Bl for a fixed set of atomic 
chemical potentials [see Sect. Ill Dl . (//) The concentrations of 
the intrinsic charge carriers and the intrinsic defects are eval- 
uated using Eqs. (O and ©■ (Hi) The conductivity is calcu- 
lated as described in Sect. Ill CI and the oxygen chemical po- 
tential is converted to an oxygen partial pressure according to 
Eq. ( TToT ). In the following we explicitly assume that the mate- 
rial is always able to reach equilibrium, which is a reasonable 
assumption at elevated temperatures. Further computational 
details are given in the appendix. 



III. RESULTS AND DISCUSSION 

A. Equilibrium conductivity at elevated temperatures 

We have implemented the model described above and used 
the formation energies from Ref. |H A band gap of Eq = 
3.0 eV was employed which is 0.4 eV smaller than the zero 
Kelvin-extrapolated band gap mimicking the shrinking of the 
band gap with increasing temperature. The value of 3.0 eV 
lies between the value obtained by temperature scaling of the 
band gap reported by Wemple,- which yields approximately 
2.8 eV at 1400 K, and the values for the band gap discussed 
by Chan et alJ^, which range between 3.0 eV and 3.4 eV. 
Both undoped and weakly ("accidentally") doped materials 
( n< A = 10 17 cm_3 ar, d 10 18 cm -3 ) were considered. The cal- 
culated equilibrium conductivity as a function of temperature, 
impurity concentration and oxygen partial pressure is shown 
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FIG. 2. (Color online) (a) Calculated conductivity as a function 
of oxygen chemical potential along the line A-D (i.e., for barium- 
rich conditions) in Fig. [TJ The results in the absence of any impu- 
rities (tia = 0cm~ 3 ) are shown by dotted lines; solid and dashed 
lines correspond to acceptor doping levels of tia = 10 17 cm -3 



and 
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10 cm , respectively, (b) Comparison of calcu- 



lated (thick solid lines) and experimentally measured conductivity 
curves (thin lines and symbols). Experimental data from Refs. \V3- 
[l5l For the calculations in (b) an accidental acceptor doping level of 



n A 



was adopted. 



in Fig. [2] All curves display the shape characteristic for a tran- 
sition from n-type (negative slope) to p-type (positive slope) 
conduction. 



1. Undoped material 

First we consider an ideally pure material (n™ 1 = cm -3 ) 
under Ba-rich conditions (i.e., for chemical potentials along 
A-D in Fig. [l} for which the thin dotted lines in Fig.f2]are ob- 
tained. Throughout the n-type region a slope of — 1/6 is ob- 
served [compare the dashed line segments in Fig. [2 a)], which 
changes to +1/6 in the p-type region. Analysis of the defect 
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(a) pure, 1 073 K, Ba-rich 
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(b) rf/ = 10 17 cm" 3 , 1073 K, Ba-rich 
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FIG. 3. (Color online) Charge carrier and defect concentrations 
for (a) pure and (b) accidentally acceptor doped material (n™ = 
10 17 cm -3 ). 



concentrations [Fig. [2 a)] shows that in the n-type region dou- 
bly charged oxygen vacancies are the dominant defects which 
gives rise to a slope of — 1/6. This slope can also be derived if 
one treats the point defect equilibria in the material using the 
mass-action law approach. Starting from the point defect 
reaction 

00 V + 2e' + ^0 2 , 

one obtains a mass-action law which links the oxygen partial 
pressure, p 02 to the concentration of doubly charged oxygen 
vacancies [V ], and the concentration of electrons n e 



Ki = [v ]n 2 J£. 



(17) 



Application of the Brouwer approximation (one defect domi- 
nates the charge neutrality condition, Eq. (01) for the oxygen 
vacancies [V ] = 2n e then yields 



-1/6 

n e ccp J . 



(18) 



which because of Eq. ( fTTT i leads to the same slope in the con- 
ductivity. 
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FIG. 4. (Color online) Variation of the electron chemical potential 
with both temperature and impurity concentration. The curves move 
downward, i.e. to more p-type conditions, with decreasing temper- 
ature as well as increasing impurity concentrations. The gray bars 
indicate the reduction of the band gap at higher temperatures. 



In the p-type region (Vn — V )" di-vacancies dominate, 
leading to a slope of +1/6. The latter can be derived using 
the mass-action law approach as follows. The point defect 
reaction 



Ba0 + ^0 2 <-> (Vr 



F )" + 2/i+BaTi03 



leads to 



"1^2-1/2 



A u = [(Vh - V )"]nip 



(19) 



which using the simplified charge neutrality condition [(Vn 
V )"] = 2n h yields 



+ 1/6 



(20) 



For clarification it should be pointed out that if the con- 
ductivity is plotted not against the oxygen partial pressure but 
against the oxygen chemical potential the minima all occur at 
the same chemical potential. The gradual shift of the minima 
in Fig. [2]is thus merely a consequence of the temperature de- 
pendence of the relation between the chemical potential and 
the oxygen partial pressure described by Eq. ( IT6l ). 

The self-consistently determined electron chemical poten- 
tial is shown in Fig. [4] At low oxygen partial pressures it is 
located in the upper half of the band gap corresponding to 
?i-type material whereas for larger oxygen partial pressures 
the electron chemical potential resides in the lower half of the 
band gap. Figure|4]shows that for an ideally pure material the 
pressure dependence of the electron chemical potential does 
not change with temperature. This observation contrasts with 
the temperature induced shift of the minimum of the conduc- 
tivity curves fFig.|2ta)]. 

The temperature and pressure dependence of the conductiv- 
ity as well as the electron chemical potential obtained under 
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Ti-rich conditions (i.e., for chemical potentials along the line 
B-C in Fig. [TJ very closely resemble the results under Ba- 
rich conditions. In the n-type region doubly charged oxygen 
vacancies are again the primary defects. In the p-type region, 
however, doubly charged barium vacancies dominate. Follow- 
ing a similar derivation as for the Vn — Vq di-vacancies (see 
e.g., Ref.O one again obtains a slope of +1/6. 

As demonstrated above the use of the mass-action law and 
the Brouwer approximation allow to deduce the slope of the 
curves for situations in which one defect dominates. Such an 
approach is, however, bound to fail in any transition region or 
in regions where several point defects have similar concentra- 
tions. This is for example the case near the crossings of the 
lines in Fig. [3] This limitation is avoided by using the full ap- 
proach outlined in Sect.|II]which furthermore does not require 
any presumptions with regard to the prevalence of any partic- 
ular defect reaction. In addition it allows to obtain the con- 
centrations of secondary point defects which although they do 
not affect the charge carrier concentrations still can impact the 
electronic properties through carrier scattering and trapping. 



2. Weakly acceptor doped material 

If under Ba-rich conditions a low concentration of accep- 
tors is present in the material, one obtains the bold solid lines 
in Fig. |2] (n™ 1 = 10 17 cm -3 ). For low oxygen partial pres- 
sures and high temperatures they have a slope of — 1 /6 in the 
?i-type region just as in the undoped material. Again this is 
due to positively charged oxygen vacancies as illustrated in 
Fig. 0b). As the temperature is lowered and/or the oxygen 
partial pressure rises, a transition to a slope of — 1/4 is ob- 
served. As shown in Fig. |5tb) this change corresponds to 
the onset of extrinsic behavior, i.e. the dominant source for 
holes are no longer intrinsic but extrinsic defects. In this case 
the simplified charge neutrality condition reads [V@] = 2n e ^- 
which if inserted into ( TPTI i also leads to a slope of +1/4. At 
higher oxygen partial pressures and higher temperatures the 
slope changes to +1/6 — as in the ideal case — indicating in- 
trinsic behavior and dominance of Vn — Vq di-vacancies. In 
contrast, at lower temperatures a slope of +1/4 is observed 
which is consistent with extrinsic behavior^ 

As shown in Fig. [4] the electron chemical potential as a 
function of the oxygen partial pressure again displays a re- 
duction of the slope for pressures > 10 -3 atm which results 
from the coupling of the concentrations of intrinsic holes and 
extrinsic defects which are of similar magnitude in this range. 
For even larger oxygen partial pressures [outside the range of 
Fig.[2]but visible in Fig.[3jb)], the concentration of Vr — V@ 
di-vacancies exceeds the concentration of extrinsic acceptors 
and the slope of the conductivity curves reverts to +1/6. The 
transition from n-type to p-type conduction correlates with a 
significant variation of the electron chemical potential over the 
band gap. In contrast to the case of an ideally pure material, 
material which contains extrinsic defects exhibits a marked 
temperature dependence of the electron chemical potential vs 
pressure curves. As the temperature is reduced the electron 
chemical potential curves are pushed downwards, which leads 



to the remarkable finding that for a certain range of oxygen 
partial pressures, the electron chemical potential moves from 
the upper to the lower half of the band gap as the temperature 
is lowered, indicating a transition fro n to p-type conduction. 

If the dopant concentration is further raised (thin dotted 
lines in Fig. [2] n™ 1 = 10 18 cm -3 ) extrinsic acceptors dom- 
inate over the entire range of chemical potentials and a slope 
of —1/4 (+1/4) is obtained throughout the n-type (p-type) 
region. The temperature dependence of the electron chemical 
potential curves (Fig. |4j is even more pronounced than in the 
case of rf^ = 10 17 cm~ 3 . 

The entire situation is very similar if Ti-rich conditions are 
imposed, the major difference being again the occurrence of 
barium vacancies instead of Vn — V@ di-vacancies in the p- 
type region. 

B. Comparison with experiment 

The conductivity of both nominally undoped as well as in- 
tentionally doped barium titanate has been repeatedly mea- 
sured as a function of oxygen partial pressure and at ele- 
vated temperaturesi&^ r 15 - 30 ' 31 These studies provide a com- 
prehensive data set for comparing our calculations with ex- 
periment. In Fig. |2jb) the results of several measurements are 
plotted together with the curves calculated for a doping level 
of rffr = 10 17 cm -3 . The agreement is very good. The cal- 
culations reproduce the n-type/p-type transition, the tempera- 
ture dependence of the position of the minima as well as the 
changes in the slopes. 

In order to explain the experimentally observed transition 
within the n-type region from a slope of —1/6 to a slope of 
— 1/4, two different models have been discussed: (/) The most 
early studies proposed the transition to be related to a change 
of the charge state of the oxygen vacancy^ The measure- 
ments could be reproduced using a model in which the oxy- 
gen vacancy +1/ + 2 transition level is located about 1.3 eV 
below the CBM and thus very close to the center of the band 
gap. (if) Most studies (see e.g., Refs. [l3i[l4ll3ll andHH), how- 
ever, assume that even in the most carefully prepared samples 
a background concentration of "accidental" acceptor impuri- 
ties is present which gives rise to the transition between the 
slopes. 

The present calculations in conjunction with the DFT data 
from Ref. HI provide very strong evidence for the second ex- 
planation. In order to obtain further support for this picture, 
we artificially pushed the +2 / + 1 transition level, which is 
located just 0.05 eV below the conduction band minimum, to- 
ward the the middle of the band gap by reducing the formation 
energy of the singly charged oxygen vacancy. The thus ob- 
tained conductivity curves do indeed display a transition from 

a oc J*0 1 ^ 6 to a oc "p^J^ ■ However, in order to reproduce at 
least approximately the experimental data the formation en- 
ergy of Vq had to be reduced by about 1 eV which is signifi- 
cantly larger than the error bar of the DFT calculations. 

The p-type region within which the slope is positive is dom- 
inated by acceptor defects. Since experimentally one observes 
a slope of about +1/6 in this region, it has been widely as- 
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FIG. 5. (Color online) "Phase diagram" illustrating the prevalent 
defects as a function of temperature and oxygen partial temperature. 
The thick solid lines separate the regions within which the indicated 
defects dominate. The dotted and dashed lines connect points along 
which the concentration of the dominant defect is constant (dash-dot- 
dot: 10 18 cm" 3 , dash-dot: 10 17 cm -3 , dotted: 10 16 cm" 3 , dashed: 
10 15 cm" 3 ). 



sumed that barium vacancies in charge state —2 are responsi- 
ble for this behavior. If one considers single vacancies only 
the alternative intrinsic acceptor defect would be the titanium 
vacancy, which occurs in charge state —4 and thus would lead 
to a slope of +1/5. The discussion above, however, shows 
that at least under Ba-rich conditions the dominant defect is 
the Vri — Vq di-vacancy, which — equivalent to the barium 
vacancy — gives rise to a slope of +1/6. Thus, on the basis 
of the conductivity curves alone the intrinsic acceptor defect 
cannot be determined unambiguously. 



C. Equilibrium defect concentrations 

The successful validation of our calculations through com- 
parison with experimental data demonstrates the capacity of 
DFT calculations and allows to use the present calculations 
for obtaining a more detailed picture of the thermodynami- 
cal behavior of point defects in this material. We can thus 
determine the prevalent intrinsic point defects depending on 
temperature and oxygen partial pressure. The results obtained 
for an ideally pure material are shown for Ba and Ti-rich con- 
ditions in Fig. [5] In this figure thick solid lines confine the 



regions within which a certain defect prevails. The concen- 
tration of the dominant defects is shown by the dashed lines 
along which the concentration is constant. 

The diagrams can be read as follows: Assume a mate- 
rial is synthesized at HOOK under Ba-rich conditions and an 
oxygen partial pressure of 10~ 2 atm [point A in Fig. |5jb)]. 
The material contains an acceptor impurity concentration of 
10 16 cm~ 3 . Under the synthesis conditions the dominant de- 
fect is the oxygen vacancy in charge state +2. The material 
is subsequently annealed at a temperature of 1200K while 
keeping the oxygen partial pressure at 10" 2 atm [point B in 
Fig. Ob)]. Under these conditions the dominant defect is the 
Vri — Vq di-vacancy. If the material is cooled below about 
1000 K the dotted line in Fig. |5jb) indicates that the concen- 
tration of — Vgj di-vacancies falls below 10 16 cm" 3 [point 
C in Fig. 13b)]. Therefore, at temperatures below 1000 K 
the accidental acceptor dopants dominate the charge equilib- 
rium and the material displays extrinsic behavior [point D in 
Fig. |3b)]. It is important to point out that this analysis is 
strictly valid only in thermodynamic equilibrium. 



IV. SUMMARY AND CONCLUSIONS 

By combining thermodynamic considerations and several 
basic relations of semiconductor physics we have obtained 
a concise scheme for the modeling of electrical conductivi- 
ties on the basis of first-principles calculations. Compared 
to "classical" defect models which are based on mass-action 
laws to connect defect concentrations, the present scheme re- 
quires a minimum number of approximations. In particular, it 
does not rely on any assumptions with regard to the prevalence 
of any particular defect reaction or defect. 

We have applied this scheme to BaTiOa, which is impor- 
tant both from the technological and the fundamental perspec- 
tive, using a complete set of thermodynamic data on intrin- 
sic point defects obtained from density-functional theory cal- 
culations. A numerical algorithm was implemented for the 
self-consistent determination of the electron chemical poten- 
tial, which enabled an extensive analysis of the dependence of 
the electrical conductivity on chemical environment as well as 
temperature. 

In agreement with earlier experimental studies our analy- 
sis has shown that the n-type conductivity, which is observed 
under low oxygen partial pressures, is due to doubly charged 
oxygen vacancies. The p-type region, which is observed at 
larger oxygen partial pressure, is caused by barium vacancies 
and Vji — V® di-vacancies under Ti and Ba-rich conditions, 
respectively. It needs to be stressed that since both of these 
defects occur in charge state —2 and therefore lead to a slope 
of +1/6 in the conductivity vs oxygen partial pressure plot, 
they cannot be distinguished on the basis of the conductivity 
curves alone. 

Our approach furthermore allows us to determine the evo- 
lution of the defect concentrations under "true" equilibrium 
conditions, i.e. in the absence of any kinetic barriers. We have 
employed this possibility to establish a point defect "phase di- 
agram" which displays the dominant intrinsic point defect as 
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a function of temperature and chemical environment. 

The successful application of the scheme outlined in the 
present paper demonstrates the predictive power of first- 
principles calculations. It also constitutes the stepping stone 
for future work which should address the effects of kinetic 
barriers and implement a more complex treatment of extrinsic 
defects. 
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Appendix A: Computational details 

1. General remarks 

The density functional theory (DFT) calculations^ from 
which we obtain our input data provide values for the ener- 
gies of formation E^. In order to obtain the free energies of 
formation AG{ the (vibrational) entropies of formation AS{ 
and the formation volumes AV/ are required 

AG{ = E( - TAS{ + P AV/. (Al) 

Since we are interested in ambient pressures (i.e., p « 0) the 
last term is virtually zero. It is in principle possible to deter- 
mine the vibrational entropy but it requires very large super- 
cells which is computationally extremely demanding ^ In the 
present work we, therefore, simply set all defect formation 
entropies to which at 1200 K amounts to a reduction of 
the free energy of formation by 0.16eV compared to the zero 
Kelvin value. 



2. Effect of band gap corrections 

The formation energies in Ref. [H were calculated within 
density functional theory using the local density approxima- 
tion. Since this calculation method is subject to a substantial 
band gap underestimation, a correction scheme was applied 
which implements a rigid shift of the conduction versus the 
valence band states. In simple terms, the difference between 
the experimental Eq P and the calculated band gap Eq 1c is dis- 
tributed between the valence and the conduction band states, 

A&g = - Ef c = AE WB +AE CB . (A2) 

Unfortunately, the ratio of A£Vb and AEqb is unknown. 
Even the GW-method^ which in principle is capable of pro- 
viding this information and which works well for many non- 
oxide materials, fails and yields a considerable overestima- 
tion of the band gap^ In Ref. [5] we therefore simply as- 
signed the band gap error entirely to the conduction band, 



AEcb = AEq~, AEy B = 0. This choice has neither an im- 
pact on the location of the equilibrium transition levels nor on 
the conclusions in Ref.d. The values of A£Vb an d AEcb do, 
however, affect the absolute values of the formation energies 
and are therefore important in the present work. 

By explicit calculation one can show that the effect of shift- 
ing A_Evb versus Ai?cB is equivalent to rigidly shifting the 
conductivity curves in Fig. [2] (details below) along the pres- 
sure axis. Neither the shape, the slopes, nor the magnitude of 
these curves are affected. We have therefore decided to adjust 
the ratio of A£Vb ar, d A£cb such that the minimum of the 
conductivity at the highest temperature considered (1473 K) 
is located at the same oxygen partial pressure as in the exper- 
iments. The final values are AE WB = OAbAE^' = 0.76, eV 
and AE CB = 0.55A^g r = 0.92 eV, which are actually of a 
very reasonable magnitude, considering that a rule of thumb 
for many semiconductors is a ratio of 2: 1 for AEqb'AEvb- A 
shift of 0.1 eV in either direction amounts to a shift in along 
the pressure axis by 10 ±2 atm. All data discussed in the fol- 
lowing was obtained using the values for A£Vb and A£cb 
quoted above. 
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